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ABSTRACT 


Two  problems  involving  the  derivation  of  bounds  on  distribu¬ 
tions  with  a  decreasing  failure  rate  (DFR  distributions)  are  presented. 
Given  that  an  item  has  a  decreasing  failure  rate,  sharp  upper  and 
lower  bounds  on  the  burn-in  time  to  achieve  a  specified  mean  residu¬ 
al  life  are  derived  The  bounds  rely  only  on  the  DFR  assump¬ 
tion  and  a  knowledge  of  the  first  moment  and  a  percentile  of  the 
failure  distribution. 

An  early  estimate  of  the  five  year  survival  proportion  (com. 
monly  called  the  five  year  cure  rate)  is  of  great  interest  in  assessing 
the  value  of  a  treatment  for  a  mortal  disease  such  as  cancer.  As¬ 
suming  that  the  distribution  01  time  to  death  is  DFR  and  assuming 
a  knowledge,  of  the  mean  and  a  percentile,  sharp  upper  and  lower 
bounds  on  the  survival  proportion  are  obtained. 

In  addition  some  bounds  on  the  hazard  rate  and  density  of  a 
DFR  distribution  are  given. 


1 


1.  Introduction.  The  need  for  complex  electronic  equipment  in  lo¬ 
cations  where  replacement  of  failed  parts  is  impossible  (e.g.  , 
ballistic  missiles,  satellites,  etc.)  has  necessitated  the  production 
of  very  reliable  components  such  as  semiconductors.  One  method 
v'hich  is  almost  universally  used  to  help  achieve  high  reliability  is  to 
pre-age  or  burn-in  the  components  to  eliminate  the  "sports"  or  early 
failures.  It  has  long  been  known  that  semiconductors,  for  example, 
exhibit  infant  mortality  yet  do  not  wear  out  (See:  Blakemore, 
Kronson,  Von  Alven,  1963;  Norris,  1963;  Von  Alven,  1962;  and  Von 
Alven,  Blakemore,  1961);  that  is,  they  exhibit  a  decreasing  failure 
rate  (are  DFR  ).  The  problem  which  now  presents  itself  is  how 
long  to  burn-in  the  component  to  achieve  a  specified  reliability.  The 
problem  has  been  answered  in  the  past  by  assuming  a  Weibull  distri¬ 
bution  of  time  to  failure  (Von  Alven,  1962;  Watson,  Wells,  1961),  for 
which  there  seems  little  statistical  validation.  However,  this  as¬ 
sumption  is  not  necessary.  It  suffices  to  assume  that  the  distribu¬ 
tion  is  DFR  and  its  first  moment  and  a  percentile  are  known  in 
order  to  obtain  sharp  bounds  on  the  residual  mean  life.  It  is  a  simple 
matter  once  the  bounds  are  known  to  determine  the  minimum  burn-ir. 
time  to  achieve  a  specified  residual  mean  life.  This  is  shown  in 
section  2.  Although  the  burn-in  problem  is  developed  in  terms  of  a 
particular  example,  burning-in  can  clearly  be  used  to  advantage  on 
any  item  which  exhibits  a  decreasing  failure  rate. 

A  problem  which  bears  some  slight  similarity  to  the  burn-in 
problem  is  that  of  estimating  the  five  year  survival  proportion  in  a 
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population  of  cancer  patients  (also  called  five  year  cure  rate 
Berkson,  Gage  1952  '•  When  a  new  treatment  is  tested  clinica  lly 
it  is  desirable  to  obtain  an  estimate,  as  early  as  possible,  of  its  ef¬ 
fectiveness.  One  objective  indication  of  this  is  the  five  year  survival 
proportion.  Berkson  and  Cage  (1952)  have  estimated  a  related  quan¬ 
tity,  the  cure  proportion,  by  assuming  in  the  interests  of  mathemati¬ 
cal  expediency  that  the  death  rate  due  to  cancer  is  a  constant.  It  is 
well  validated  that  the  death  rate  from  cancer  is  decreasing  with 
time  (Cutler,  Axtell,  1963,  Berkson,  Gage,  1952)  and  this  fact  is 
used  together  with  a  knowledge  of  the  mean  and  a  percentile  of  the 
distribution  of  time  to  death,  to  obtain  sharp  upper  and  lower  bounds 
on  the  five  year  survival  proportion.  This  is  discussed  in  section  3. 

The  most  appealing  feature  of  the  solutions  presented  to  the 
preceding  two  problems  is  that  no  assumption  of  a  parametric  ex¬ 
pression  for  the  probability  distribution  is  made.  Such  an  assump¬ 
tion  would  be  very  difficult  to  verify  using  the  truncated  data  from 
semiconductor  .ife  tests  or  the  small  sample  data  resulting  from  a 
clinical  trial.  Yet  it  is  possible  to  obtain  a  reasonable  estimate  of 
the  mean  life  and  early  percentile  for  a  truncated  life  test.  This 
information  together  with  the  DFR  assumption  enables  us  to  obtain 
bounds  on  the  relevant  quantities  to  be  evaluated. 

Nearly  all  the  work  completed  to  date  on  bounds  for  distribu¬ 
tions  possessing  a  monotone  failure  rate  has  been  done  by  two 
authors,  Barlow  and  Marshall  (see:  Barlow,  1963,  Barlow  and 
Marshall,  1963,  !964a  and  1964b),  who  have  often  collaborated. 
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Barlow  and  Proschan  (1964)  also  develop  many  applications  of  these 
bounds  in  the  field  of  reliability  theory. 


Only  the  relevant  bounds  are  set  forth  in  sections  2  and  3t 
while  the  bounds  are  derived  in  section  4.  A  few  additional  theorems 
are  added  to  complete  the  discussion  of  bounds  for  DFR  distribu¬ 
tions. 

Mathematical  Preliminaries.  If  the  failure  distribution  F  has 
density  f  ,  then  the  failure  rate  q(t)  is  defined  for  those  values 
of  t  for  which  F(t)  <1  by 


F(t) 


where  F(t)  =  1  -  F(t)  and  it  is  assumed  that  F(0")  =  0  .  The 
failure  rate  is  also  known  as  the  hazard  rate  and  by  actuaries  as  the 
"force  of  mortality". 

It  can  be  readily  verified  that  q(t)  =  -  log  F  (t)  ,  when  a 
density  exists.  Hence  F  is  DFR  (IFR)  if  log  F  (t)  is  convex 
on  [  0,or]  (concave  where  finite).  This  fact  forms  the  basis  of 
many  of  the  proofs  of  section  4. 
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2.  The  Burn-in  Problem.  The  only  items  which  are  burnt-in  to  any 
extent  are  semiconductors,  although  burn-in  could  well  be  applied 
more  widely.  Hence,  we  use  semiconductors  to  illustrate  our  dis¬ 
cussion  of  the  burn-in  problem. 

A  few  remarks  will  made  on  the  life  distribution  of  a  semicon¬ 
ductor  since  there  is  still  some  disagreement  as  to  the  actual  shape 
of  the  distribution.  Peck  (Von  Alven,  1962,  Chapter  2)  asserts  that 
his  data  on  semiconductors  exhibits  first  an  increasing  failure  rate 
and  then  a  decreasing  failure  rate.  He  supports  this  by  a  physical 
explanation  of  the  cause  of  the  initial  IFR  .  Other  authors  (i.e., 
Norris,  1963}  maintain  the  distribution  is  DFR  and  the  IFR  ,  but 
this  theory  is  more  intuitive  than  well  supported  by  life  test  data. 

The  great  bulk  of  work  in  the  life  testing  of  semiconductors  supports 
the  theory  that  there  is  no  wear-out,  and  the  failure  rate  is  always 
decreasing.  In  particular,  ARINC  Research  Corporation  (Blake- 
more,  Kronsor.,  Von  Alven,  1963;  Von  Alven,  Blakemore,  1961) 
tested  10,  300  individual  devices  with  approximately  200,000  separate 
life  test  measurements  and  reported  they  could  not  detect  any  wear- 
out,  but  found  in  almost  every  case  the  life  distribution  was  DFR  . 

The  time  to  burn-in  to  achieve  a  specified  reliability  has  been 
determined  in  practice  by  assuming  the  life  distribution  is  either 
Weibull  or  lognormal.  This  affords  a  very  rapid  and  simple  method 
for  evaluating  the  burn-in  time  but  in  the  n^ajority  of  cases  where 
these  distributions  are  assumed,  it  is  done  with  little  statistical  val¬ 


idation. 
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\  non-paramftric  approach  based  purely  on  the  DFR  assump¬ 
tion  obviates  this  uncertainty  of  distribution  validity  and  gives  shirp 
bounds  which  although  more  conservative,  Hoes  guarantee  achieving 
the  specified  reliability  since  it  is  valid  for  a  larger  class  of  distri¬ 
butions . 

Bounds  are  set  out  below  on  the  survival  probability  and  the 
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T  he  proofs  of (2 . 1) thr ough(2 . 4) are  given  in  theorems  (4./.)  through 
(4.  h)  of  Section  4. 

The  most  interesting  result,  the  lower  hound  on  the  mean  re¬ 
sidual  life,  is  useful  only  up  to  time  £  ,  for  after  this  time  the 

hound  is  constant  and  so  gives  no  indication  of  the  effect  of  increased 
burn-in.  However,  even  in  the  light  of  this,  the  restriction  on  the 
percentile  that  <  1  or  1  -  p  <  e  P  does  not  appear  restric¬ 
tive  in  the  burn-in  problem  as  it  is  unlikely  that  the  item  would  be 
burnt-in  for  a  time  greater  than  its  mean  life. 

By  the  use  of  equation  (2.  3),  the  burn-in  time  to  achieve  a 
specified  mean  residual  life  has  been  calculated  and  is  given  in 
tables  1  through  3. 

Bounds  in  terms  of  the  first  and  second  moment  have  little 
usefulness  since  the  sharp  lower  bound  on  th'*  mean  residual  life  in 
this  case  is  the  value  of  the  mean,  thu'j  giving  no  indication  of  the 
benefit  gained  by  burn-in.  Bounds  on  the  survival  probability  for 
this  case  have  been  calculated  by  Barlow  and  Marshall  (1963). 

The  burn-in  of  semiconductors  is  sometimes  carried  out  at  an 
increased  stress  to  accelerate  the  ageing,  (free  Norris,  1963).  This 
may  introduce  many  complications  in  the  prediction  of  the  optimal 
burn-in  time.  However,  it  will  be  assumed  that  it  is  possible  to  con¬ 
vert  from  the  un-accelerated  burn-in  time  to  the  accelerated  burn-in 


time.  This  could  possibly  be  done  with  the  use  of  the  Arrhenius 
relation,  although  v0n  Alven*  states  that  the  Eyring  rel-itionsnip 
may  possibly  be  better. 
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3.  Bounding  the  Survival  Proportion.  A  frequently  used  objective 
index  of  the  effectiveness  of  a  treatment  is  the  proportion  of  patients 
surviving  the  disease  for  five  years.  (See  for  example,  Berkson, 
Gage,  1952).  This  will  be  called  the  five  year  survival  proportion 
although  it  has  been  generally  called  by  the  misleading  name  of  "five 
year  cure  rate"  by  the  medical  profession.  Clearly  the  word  cure  is 
inappropriate  as  even  in  so  mortal  a  disease  as  cancer  it  is  not  cer¬ 
tain  that  the  patient  is  cured  when  he  has  survived  five  years. 

Berkson  and  Gage  (1952),  have  discarded  the  idea  of  a  survival 
percentage  in  favour  of  computing  the  cure  proportion;  defining  the 
patient  cured  when  his  death  rate  is  the  same  as  the  normal  mortali¬ 
ty  rate.  Catler  and  Axtell  (1963)  point  out  that  in  some  cancers  this 
is  never  achieved  and  they  redefine  cure  as  the  achievement  of  a 
stabilized  death  rate.  But  with  both  Berkson,  et  al  and  Cutler, 
et  al,  information  is  required  over  a  long  span  of  time  to  estimate 
the  cure  proportion. 

Although  the  five  year  survival  proportion  does  not  give  as  ac¬ 
curate  a  picture  as  the  cure  proportion  (where  the  latter  is  relevant) 
it  does  afford  a  good  indication  of  the  effectiveness  of  a  treatment 
and  can  be  calculated  more  simply  and  at  an  earlier  time  than  can 
the  cure  proportion. 

An  extension  of  the  work  already  done  in  this  field  would  be  to 
predict  the  five  year  survival  proportion  at  the  end  of  only  say  one 
year  of  a  -.niral  tria.  thus  enabling  an  ear’y  assessment  of  the 


treatment  to  be  made. 
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The  mathematical  model  will  be  simplified  by  assuming  that 
the  probability  of  death  due  to  normal  causes  is  independent  of  the 
probability  of  death  due  to  cancer.  This  clearly  oversimplifies  the 
issue  but  Berkson  and  Gage  (1952'  maintain  that  this  assumption  does 
give  reasonable  results.  Thus 
(3.1)  F(t;  ”  Fn(OF(t) 

where  F  (t)  is  the  survival  probability  in  time  t  and  F^(t)  and 

FJt)  are  respectively  the  probabilities  of  death  from  "normal 

causes"  and  from  cancer  in  time  t  .  It  is  assumed  that  F  is 

n 


known  from  life  tables.  Also 
(3.2)  q(t)  =  qc(t)  +  qn(t) 

where  q(t)  is  the  death  rate  at  time  t  and  the  subscripts  c 
and  n  are  as  in  (3.1). 

Berkson  and  Gage  (1952)  have  shown  that  qc(t)  is  decreasing 
at  a  rate  which  is  a  function  of  the  mean  time  to  death  of  the  untreat¬ 
ed  patients.  The  normal  mortality  rate,  q^(t)  ,  is  increasing. 

Thus  q(t)  is  initially  decreasing  and  then  increasing. 

If  the  time  at  which  q(t)  changes  from  decreasing  to  in¬ 
creasing  is  large  compared  with  5  years,  bounds  on  Fc  (5  years) 
may  be  estimated  by  assuming  that  F  is  DFR  .  Thus  from  a 
knowledge  of  the  mean  and  a  percentile  of  F  ,  bounds  on  F  (5 
years)  may  be  determined  by  using  equations  (2.1)  and  (2.2),  and 
from  this  F^  (5  years)  can  easily  be  deter  mined  fromequation  (3. 1). 

A  less  accurate  but  surer  method  would  be  to  adjust  the  ob¬ 
servations  of  the  time  to  death  to  obtain  an  estimate  of  the  DFR 
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distribution  F,  .  Hence  by  a  straightforward  application  of  the 
DFR  bounds,  we  may  obtain  bounds  on  (5  years). 

Bounds  on  the  failure  rate  at  five  years  may  also  be  evaluated 
to  give  a  further  indication  of  the  effectiveness  of  the  treatment. 
These  bounds  can  be  estimated  by  the  use  of  the  DFR  bounds 
given  below  and  equation  (3.  2) 

F  DFR  ,  -  I  ,  F(Cp) ’  =  p  and  5p<  1  or  1  -  p<  P 

(3.3)  q(t+)  >  a 

where  a  is  defined  by  the  unique  solution  to  (2.5). 

r  00  t  =  0 

(3.4)  q(t")  <  -t 

L  b1  t  >  0 

See  theorems  (4.6)  and  (4.7)  for  proofs. 
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4.  Bounds  on  DFR  Distributions. 

4.1.  Bounds  Given  the  Mean  and  a  Percentile  have  great  ap¬ 
peal  since  good  estimates  of  these  two  quantities  can  be  obtained 
from  limited  life  test  data.  In  addition,  the  desired  bounds  are  all 
nontrivial  which  is  in  contrast  to  the  case  where  the  first  two  mo¬ 
ments  pj  and  p^  are  given.  In  this  case  the  lower  bound  on  the 
mean  residual  life  is  trivial  and  sharp  (i.e.  is  the  value  p^)  . 

The  DFR  distribution  with  mean  p^  =  1  implies  an  upper 
bound  on  the  percentile.  The  lower  bound  is  zero,  and  both  bounds 
are  sharp. 

r  e'?p  ,  €  <  1 

(4  1)  0  <  1  -  p  <  J  P 

1  <*«>  >  1 

P  .  P 


This  result  is  proved  by  Barlow  and  Marshall  (1963). 


Two  distributions  will  be  defined.  Let 


(4.2) 

f0e‘blX  .  x  <  ? 

7b/x*  =  L.  .  -(x-$  )b2  . 

1  4l-p)e  '  'p'  2  ,  x  >  €p 

where  7,  = 

bi 

1 

cr 
►— • 

and  /3  and 

b^  are  given  by 

(4.3) 

0  e”bl?P  =  1  -  p  ,  (0  <  0  <  1) 

and 

C  X  (x)  dx  =  1  ; 

Jo  °1 

Clearly  J, 
bi 

has  mean  p^  =  1  and  ptb  percentile 

be  the  subclass  of 


1 


Let 


which  is  DFR  .  A 
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lemma  will  now  be  proved  to  facilitate  the  proof  of  theorem  4. 1. 


Lemma  4.1  If  the  percentile  £  of 

-  K  p 

ditions  (4^1)  then  the  distribution  of  the  class 
maximum  value  to  the  mass  at  the  origin  is 


\ 

\ 


satisfies  the  con- 

which  gives  the 

* 

,  where  b^  is 


the  minimum  value  of  b^  such  that  b^  =  b^  ;  i.e.,  an  exponential 
on  [0,<r]  ,  with  possible  mass  at  the  origin. 


Proof.  Clearly  the  distribution  of  the  class  which  attains 

maximum  mass  at  the  origin  is  given  by  b^  =  0  ;  but  it  is  not  DFR 
since  it  must  be  log  concave.  As  b^  increases  from  zero  (thus 
decreasing  /3  )  ,  b^,  increases  but  less  rapidly,  and  the  distribu¬ 
tion  will  first  be  DFR  when  b^  =  b^  ,  attaining  then  the  maxi¬ 
mum  mass  at  the  origin  for  a  distribution  R,  .  The  equation  of 

1 

this  distribution  is 

(4.4)  \=-b2<X)  =ae_*X 


where 

(4.5) 


a  is  the  minimum  solution  of 
1  -  p  =  ae 
a  <  1  J 


Note  that  a  solution  a  is  guaranteed  by  the  conditions  on  the  per¬ 
centile,  and  that  for  ?  <  1  the  solution  to  (4.  5)  is  unique,  but  for 

£  >1  there  may  be  two  solutions.  || 


Theorem  4.1.  If  F  is  DFR  with  mean  =  1  ,  and  p*^  per¬ 
centile  £  is  given  and  satisfies  the  condition  1  -  p  <  e  ,  then 
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F  (x)  >  oe~OX  , 

x  <  « 

F (x)  <  at  OX  , 

x  >  ? 

and  a  is  defined  as  the  unique  solution  of  ere  °^P 


1  -  p  . 


Proof.  Suppose  that  the  theorem  is  false.  Then  by  convexity,  either 
(i)  or  (ii)  is  true. 


(i) 

F  (x)  <  oe'OX 

,  for  some 

X 

A 

vn 

T) 

(ii) 

F  (x)  >  ae*OX 

,  for  some 

x  >  Cp  and  F (x)  >  o-e  °X 

for 

x-  5p  • 

Suppose  case  (i): 


Figure  4. 1 


By  log  convexity  F (x)  >  a-e  aX  when  x  >  f  .  Construct  two 
exponential  curves,  e  and  e  where  f^(x)  and 

f^x)  are  linear  in  x  and  such  that 
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Now  clearly  e  >  0e  ^  for  x  >  ?  and  hence 

-Ctx 

e  1  <  Qe  for  x  <  £  .  Thus  a  DFR  distribution  has  been 

P 

constructed  with  the  same  mean  and  percentile  as  the  exponential 
K(x)  =  are  .  but  a  greater  mass  at  the  origin,  which  by  lemma 
4.  1  is  impossible. 

Suppose  case  (ii): 


Figure  4.  Z 

By  log  convexity  of  DFR  distributions  there  must  be  a  single 
crossing,  say  at  x^  ,  where  x^  >  5  .  By  assuming  a  percen¬ 

tile  at  Xj  ,  and  noting  that  the  solution  a  is  unique  since 
1  -  p  <  e  ,  it  can  be  seen  from  case  (i)  that  case  (ii)  is  im¬ 
possible.  | 

Note  that  the  condition  on  the  percentile  1  -  p  <  e  P 
which  it  was  necessary  to  assume  to  assure  a  unique  value  of  o  is 
always  satisfied  by  a  DFR  distribution  for  £  <  1  . 
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It  does  not  appear  likely  that  this  restriction  will  limit  the  applica¬ 
bility  of  the  bounds  as  it  is  hard  to  imagine  estimating  a  percentile 

at  £  >1  in  a  distribution.  In  order  to  prove  some  theorems  on 

P 

bounds,  a  class  of  DFR  distributions  Cj  is  posited  with  mean 
Pj  =  1  and  ptb  percentile  satisfying  1  -  p  <  e  . 

T  < 

(4.6) 


!  ~blX 


Ct(X)  =  le-b2x  +  <b2-bl>T 


0  <  x  <  T 
x  >  T 


where  bj  and  b^  satisfy 


(4.7) 


,.e-b,T  e-blT 


(4.8) 


1  -  P  =  exP(-b2?p  +  (b2  -  tyT ) 


and  for  this  class  to  be  DFR  ,  b^  <  b^  . 


T  >  € 


(4.9) 


-bix 

/  C  1 

CT(x)  =  {  -bzx  +  (b2-b!)T 


,  0  <  X  <  T 

,  x  >  T 


b^  and  b^  satisfy  (4.7)  and 


(4.10) 


1  -  p 


=  e‘bl5l 


and  b  «.  b  for  to  be  DFR  . 


Lemma  i  i  Assuming  1  -  p  <  e  SP  for  every  T  <  £p  there 
is  a  solution  of  (4.7)  and  (4.8)  continuous  in  T  ,  and  for  every 
T  >  5  there  is  a  solution  of  (4.7)  and  (4.10)  continuous  in  T  . 


Proof:  T  <  Cp  :  By  substituting  for  b^  from  (4.7)  into  (4.8)  we 
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obtain 


1  -  p 


exp 


{-blT-e-blT(5p-T)(l+^^i)-1| 


and  it  is  desired  to  show  this  always  has  a  solution  .  Since  the 

distribution  is  IFR  and  1  -  p  <  e  ,  we  see  b^  >  1  . 

Let 


h(b{) 


Now 


"  b  T 

expj-^T-  e~blT( £p  -  T)(l  +  e  ^  1  )'1  j  -  1  +  p 

lim  h(b.)  =  e^P-l  +  p>0 

v1  1 

lim  h(b.)  =  -1  +  p  <  0 

l_ ^ 


Thus  as  the  function  is  clearly  continuous  in  b^  ,  there  exists 
a  solution  such  that  h(b^)  =  0  . 

Since  bj  >  1  and  the  mean  =  1  ,  must  cross  e  x 

once  from  below  and  thus  b^  <  b^  . 

T  >  £  :  From  (4.10)  it  can  be  seen  that  there  is  always  a 

_ E 

solution  for  b^  . 

i.  e.  b1  =  -  »  log  (1-p)  . 

*P 

1  -  exp  f  ^  log  (1  -  p) 

Hence  b,  =  exp  ^  log  (1  -  pn  / 1  + - • - 2 - 

V*P  /l  Cp  log  (1-P) 

which  has  a  solution  continuous  in  T  for  all  T  >  ? 

P 

Since  1  -  p  <  e  ^P  and  =  1  it  can  clearly  be  seen  that 
b2  1  bl  -1 

Note  that  lim  G_(x)  =  o-e  °X  where  a  is  defined  uniquely 
T-^0  T 

by  (4.5).  This  limit  will  be  denoted  by  Gq(x)  • 
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Use  is  made  in  the  following  theorems  of  the  fact  that 
and  F  must  cross  at  least  once  if  they  have  the  same  first  mo¬ 
ment.  The  various  possibilities  for  the  intertwining  of  log 
and  log  F  are  shown  below. 
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and  the  inequalities  are  sharp,  is  defined  by  equation  (4.10) 

and  a  by  equation  (4.5). 


Proof: 


0  <  t  <  4p  From  the  log  convexity  of  F (x) 


F(t)  <  G-  (t)  . 

t  >  £  From  theorem  4.J., 

_ P 

F(t)<G0(t)  . 


0  <  t  <  £ 

-  -  P 


t  >  5. 


th 


Theorem  4.  3:  If  F  is  DFR  and  has  ^  =  1  and  p  percentile 


F  (t)  > 


r  oe 


5  given  and  also  1  -  p  <  e  ^P  ;  then 
P 

-  at 


e 


-bit 


0  <  t  <  £ 

_  _  p 

t  >  £ 


where  b^  is  defired  by  (4.10)  and  a  by  (4.  5).  The  bound  is  sharp. 


Proof: 


0  <  t  <  £  This  bound  follows  directly  from  theorem  4.1. 

~  ~  P 


£  <  t  <  oc  : 

P  _ 


The  proof  is  obvious.  The  bound  is  "epsilon" 


attained  by  lim  G_(x)  . 

T—  oc  1 


th 


Theorem  4.4:  If  F  is  DFR  ,  pj  =  1  and  the  p  percentile  £ 

-F 

given  and  also  1  -  p  <  e  P  ;  then 


f. 


F  (x)  dx 


1  .  bjt  1 
-T-  )  e  1  +  r- 


F  (t) 


>  i 


*->?p 


II 


where  is  defined  by  equation  (4.10)  and  b^  is  defined  by 

the  equations  (4.7)  and  (4.10)  with  T  =  . 


Proof: 


For  t  <  §  Since  F  (x)  <  G-  (x)  .  x  <  £ 

_ ~  P  -  *p  -  P 


r  F(x)dx  >  f 


G  (x)  dx  ,  t  <  £ 


^t  ^p 


-  P 


therefore 


^  F  (x)  dx  (  G_  (x)  dx 

Jt  >  "t  *p 


F(t) 


C.  (t) 
P 


/I  1  v  bit  A  1 

=  “"E,**  1  +B, 


For  t  >  £  The  proof  will  be  treated  in  two  cases  based 
P 

on  the  crossings  of  G_  (x)  and  F  (x)  for  x  >  £ 

Cp  -  P 


Case  (i): 


Consider  the  extremal  distribution  G-  (x)  and  let 

*P 


it  cross  F(x)  at  u(£^)  >  €  •  Obviously  due  to  log  convexity  the 

crossing  is  from  above. 


Sine  e  F  l  x)  <  G^_  (x) 

P 

n<r-  ,  .00 


I 


F  (x)  dx  i  Gr  (x)  dx 

^  Jt  ?p 


F(t) 


C,  (t) 


x  <  u (€  ) 
P 


t  <  U(?p) 


Now  by  Theorem  4.1  F(x)  <  Gn(x)  for  x  >  £  .  Thus  due 

-  U  -  p 

♦o  log  convexity  of  F(x)  ,  for  every  x  ,  u(  £  )  <  x  <  °c  ,  it  is 
always  possible  to  find  a  value  T  ,  0  <  T  <  £^ 
u(T)  =  x  . 


,  such  that 


L  3 


£V(x)dx 

F(t) 


>  inf 

T  <  € 

-  P 


f  CT(x)dx 


CT(t) 


’  tlu{?p) 


Case  (ii): 


be  shown 


u(^)  =  Cp  •  By  a  similar  reasoning  to  case  (i)  it  may 


.or 


F (x) dx  \  CT(x)  dx 

t  .  .  ,  Jt  ,  t  >  C 

-  >  inf  -  —  P 


F(t)  T<4  G_,(t) 

-  p  T 


-  inf  }■ 

T  <h 
-  P 


1  * 

=  — .  where  b,  is  the  value  of  b, 

,  *  c  L 

hl 


calculated  at  T  =  £ 


Thus,  the  lower  bound  for  the  mean  residual  life  while  not 
trivial  for  t  >  is  nevertheless  not  dependent  on  t  and  so  can 
give  no  information  on  the  merit  of  continued  burn-in. 


t  h 

Theorem  4.5:  If  F  is  DFR  and  =  1  and  the  p  percentile 

-£ 

€  is  given  and  also  1  -  p  <  e  P  ;  then 


F(x)d: 


)  dx 


F(t) 


G0(t) 


l_ 

a 


t  <  f 
-  P 


Proof:  The  proof  follows  directly  from  theorem  4.  3. 
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Theorem  4.6:  If  F  is  DFR  with  first  moment  =  1  and 

percentile  ,  and  1  -  p  <  e  ^P  ,  and  has  hazard  rate  q 

then  .  q(tf)  >  a  ,  where  a  is  the  solution  of  (4.5). 

The  bound  is  sharp. 


Proof:  The  proof  follows  directly  from  Theorem  4. 1. 


Theorem  4.7:  If  F  is  DFR  with  p^  =  1  and  percentile 

£  and  1  -  p  <  e  and  hazard  rate  q  then 

P 

f  QC  ,  t  =  0 

q(t_ )  <  i 

1  b.  t  >  0 

*  9 


where 
and  bj 


b^  is  the  solution  of  (4. 7)  and  (4.  8)  for 

is  the  solution  of  (4.10)  for  t  >  € 

-  P 


T  =  t  ,  t  <  5p  ; 
The  bound  is  sharp. 


Proof:  t  =  0  The  bound  is  trivial  and  is  attained  by  the  distribu¬ 
tion  Gq  . 

0  <  t  <  :  By  the  log  convexity  of  the  DFR  distribution  it 

can  be  seen  that  at  t 

q(t')  <  b1 

t  >  Clearly  by  log  convexity 

q(t~ )  <  bL  . 

The  bound  is  relatively  trivial  and  is  "epsilon  attained"  by  the 

distribution  lim  .  II 

T  —  oo 


.2.  Bounds  Assuming  the  First  Two  Moments  Are  Given. 


Barlow  and  Marshall  (1964  )  have  given  bounds  on  the  survival 
probability  of  the  DFR  distribution  when  the  first  two  moments  are 
given  The  same  methods  used  to  bound  the  survival  probability 


may  be  used  to  obtain  *harp  bounds  on  the  residual  mean  life,  the 
failure  rate  and  the  density  of  the  distribution.  These  bounds  are 
stated  below  and  proofs  for  them  may  be  found  in  Lawrence  (1964) 

f  F(x)dx 


(4.11) 


(4.12) 


(4.13) 

(4.14) 


no 

f  F(x)dx 
«'t 


no 


q(t')  <  a1 
q(t  )  >  az 


-  **1 


<K 

aa 


sup  a-  exp(  -a ,  t) 
T  >  t  1 


,  0<t<  1 


(415)  f(t’)  < 


max 


a  exp( -a  t) 


-  <  a  <  1 


sup 
T  >  t 


exp(-ajt) 


> .  1  <  t  < 


sup 
T  <  t 


a?exp{ -a^tf (a^-a^T } 


r  sup 


a  exp(-ot) 


-  <  a  <  1 

^2  "  - 

sup 
T  <  t 


a2exp{  -a^t  +  (a2  "a  j)T } 
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min 


inf 

-  <  o  <  1 
**2"  “ 


a  exp(-ot) 


(4  16)  f(t  +  )  >  < 


inf  a  ?exp{ -a t  +  (a^ -aj)T } 


M  2 

!>•  0  1 ' 1  7 


T  <  t 


inf  a,exp{ -a,t  +  (a  , -a. )  T) 
T  <  t 


l>7 


Where  aj 

'  a2 

are  the  unique  solutions  of  (4. 17)  and  (4.  18)  with 

T  =  t  . 

(4.17) 

1 

=  a"1  (1  -  exp(-a^T)  )  +  a^1  exp(-aJT) 

(4.18) 

^2 

=  a"2  {1  -  (a^T  +  l)exp(-aj  T)}  +  a’2  (a^T  +  l)exp( -a^  T) 

> 
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